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Abstract 



Spin-wave instability for parallel pumping is studied theoretically. The spin-wave 
instability threshold is calculated in ferromagnetic thin films under oblique field 
which has an oblique angle to the film plane. The butterfly curve of the threshold 
usually has a cusp at a certain value of the static external field. While the static field 
value of the cusp point varies as the oblique angle changes, the general properties 
of the butterfly curve show little noteworthy change. For very thin films, however, 
multiple cusps of the butterfly curve can appear due to different standing spin- wave 
modes, which indicate a novel feature for thin films under oblique field. 
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1 Introduction 

The first experiment for the parametric excitation of spin wave by parallel 
pumping was given by Schloman et al. [1]. In a parallel pumping experiment, 



a microwave magnetic field is applied parallel to the external static field. When 
the microwave field amplitude h exceeds a certain spin-wave instability ampli- 
tude h CT it, the parametric spin-wave excitation occurs. The excited spin waves 
have half the pumping frequency. The h CT i t curve plotted against the static field 
is called "butterfly curve". The spin-wave instability for parallel pumping is 
often investigated for the purpose to study relaxation phenomena in ferromag- 
netic materials, since the instability threshold depends on the balance between 
the driving power and the damping of spin waves. 
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Theoretical studies have not been very successful to explain the instability 
threshold for parallel pumping: a combination of more than two equations 
(theories) is needed to obtain the butterfly curve. In other words, even a qual- 
itative explanation for the butterfly curve cannot be successful so long as one 
tries to use only the Landau-Lifshitz (LL) equation governing the magnetiza- 
tion dynamics. By contrast the Suhl instability for perpendicular pumping has 
been well explained theoretically [2]. The Suhl instability can be described by 
using only the LL equation. For parallel pumping, however, the LL equation 
is not enough to explain the instability, and the spin-wave line width AH), is 
necessary to describe the relaxation of spin-waves as well as the instability. 

A typical butterfly curve for parallel pumping has a cusp at a certain external 
static field, which can be explained if the spin-wave line width AH^ is given. 
For example, Patton et al. proposed some trial AH^ functions and successfully 
explained butterfly curves [3] . Such a typical feature of butterfly curves usually 
does not depend on the shape of samples. For very thin films, however, more 
interesting features can appear. In fact, butterfly curves with multiple cusps 
were observed for 0.5-/zm yttrium iron garnet (YIG) films under in-plane ex- 
ternal field [4]. Those cusps were due to quantized standing- wave modes across 
the film. It is highly desirable to develop the theory of multi-cusp butterfly 
curves in the case of a general oblique field. 

In this paper, we examine the instability threshold for parallel pumping in fer- 
romagnetic thin films theoretically. The butterfly curve of the threshold will 
be calculated in cases where the external field has an oblique angle to the film 
plane. For the calculation of the threshold, we use a trial AH^ function pro- 
posed by Patton et. al. [3], which was proposed originally for spherical samples 
but can be applied to thin films [5,6]. For usual thin films, the butterfly curve 
of the threshold has a single cusp at a certain value of the static external field. 
We will show how the cusp shifts when the oblique angle changes. Moreover, 
we will reveal that multi-cusp butterfly curves can be seen for very thin films. 



2 Equations of motion for the spin- wave amplitudes 



The dynamics of magnetization field M(r) is governed by the Landau-Lifshitz 
(LL) equation, 

-d t M = MxH cS - A M x (M x H eS ). (1) 

7 M 



Here, 7 is the gyromagnetic ratio (7 < for spins); M is the value of the 
magnetization in thermal equilibrium; A is a damping constant parameter 



(A < in this case); H e g is the effective magnetic field: 

H eS = DV 2 M + H d + H a + H + h cos cot. (2) 



The first term comes from the exchange interaction; the third term is the 
anisotropy field, which is omitted below for convenience; the fourth and fifth 
terms are the external static and pumping fields, respectively. For parallel 
pumping, Hq/ /h and they are parallel to the z axis. The second term on 
right hand side of Eq. (2) is the demagnetizing filed given by the gradient of 
the magneto-static potential <ft: 

H d = -V0. (3) 



The magneto-static potential obeys the Poisson equation: 

. 47rV • M , inside the sample, 
V 2 0=<| (4) 

0, outside the sample. 



The number of independent components of M is two, since the length of the 
magnetization vector is invariant (\M | = M ). The normalized magnetization 
S = M/Mq can be represented by a point on a unit sphere. It is convenient to 
project the unit sphere stereographically onto a complex variable ip(r,t) [7]: 

D X + lOy . . 

* = TTsP (5) 

where 

^, + r q Kr - i>) q i-^* , fi x 

In terms of ip, the LL equation (1) is rewritten as 
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Fig. 1. Schematic picture of sample setting. The sample extends on the x-y plane. 
The z axis corresponds to the direction of the external filed and has an angle a to 
the z axis. 
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Now let us consider the boundary conditions, which affect the demagnetizing 
field. We assume a film of thickness d infinitely extended in the x-y plane 
as shown in Fig. 1. The external field is at an angle a to the z axis and 
corresponds to the z axis. We assume unpinned surface spins, which satisfy 
Neumann-like boundary conditions, 



9zS\z=±d/2 — 0, 



(9) 



namely, 

■% 

dz 



z=±d/2 



dz 



(10) 



z=±d/2 



Before proceeding to the linear instability of spin waves, we introduce the 
following dimensionless time and space units [8]: 



t 
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Then we obtain the linearized equations of motion corresponding to Eqs. (7) 
and (8): 



d t ij + i(l - iA) 



1 



/ 2 V 2 V - 7;( d x + i#i/)$ + (d z ® -wh- u h cosu p t)ip 



0,(12) 



V 2 $ = (d x - id y )i> + (d x + id y )ip* when -\ < z < \, 



(13) 



where 
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Now we consider the undriven case (i.e., u^ = 0), and expand ip(r,t) and 
ip*(r,t) so that they fulfills the boundary conditions (10). For even modes, 
k z = 2mn (m: integer), 



^(f , t) = J2 a fcW e^"****® cos k z z 

k 

^ (f, i) = 2 a^(t) e i (*- i+ *»») cos M. 

k 

For odd modes, fc 2 = (2m + 1)tt (m: integer) 
i]){r, t) = Y^ a k(f) e l{kxS:+ ' kv ^ sin k z z 

k 

$*(?, t) = -J2 a*_ k (t) e^+M) s i n ~k z z. 
k 



(15) 



(16) 



Here, k x = k x cos a + k z sin a, k y = k y , and /c 2 = k z cos a — k x sin a. Using the 
expansions (15) and (16), we obtain solutions of Eq. (13): for even modes, 



*(*) = -E 



r,i(fca:£+fcy2/) 



k 2 



i{k-a~ k + k + a*_~ k ) [cos fc z 5 - (-l) m e fex/2 cosh k ± z\ 

+ k z sin a ■ (a k + a^) {sin fc 2 * - (-l) m e~^/ 2 (A; z /A; ± ) sinh ~k L z}~\ , (17) 
for odd modes, 



e i(k x x+k y y) 



$(r) = - ^ ^ [K~ k - a k - MI*) {sin JM - (-l) m e~^/ 2 sinh A; ± 5J 
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K 

— k z sin a ■ (at, — a 



*_- k ) {cos fc*5 + (-l) m e-~ k ^ 2 (k z /~k ± ) cosh fc ± 5J 



where fc + = k x cosa + ik y , fc_ = fc s cos a — ifc y and /cj_ = Jk 2 + /c^. The detailed 
derivation of these solutions is shown in Appendix A. Here we approximate 
<9 2 $ = —H^/47rM with the value of uniform magnetization, k = 0: H% = 
H~ sin a + H~ cos a = —4ti(N x M x sin a + N Z M Z cos a), where N x and N z are 



demagnetizing factors. In this case, N$ = Ny = and N z = 1 • Assuming 
M z ~ M , namely, M z ~ M cosa, we obtain <9 2 $ = — H z /A-kM = cos 2 a. 
Then Eq. (12) for Uh = is rewritten as 



W + i(l - iA) 



I V ip (cosaSj + idy — sinac^)^ + (cos a — ujh)^ 



0.(19) 



Combining Eqs. (15)-(19), we obtain 

d t a k + i(l - \\)A k a k + i(l - \X)B k a*_ k = 0, 



(20) 



where 



A k = C k + ReC k , B k = (-l) n (D k e 2i ^ +C k ) . 



(21) 



Here, n is an integer, cos^fc = k x /k±, sin<y2 fc = k y /k±, k± = Jk 2 + k 2 and 



C k = cos 2 a — l 2 k 2 — ujh fk sin 2 @k 



k z 
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C fe = -^| sin a (fcj_ e 1¥,fe +k z sin a J / fe , 



Ac = -^/ fc sin 2 6» fc , 
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Sint7fc = , K z = 7TO. 
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The detailed derivation is given in Appendix B. 
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(22) 



Equation (20) represents two coupled harmonic oscillators a k , a*_ k . We diag- 
onalize Eq. (20) by means of the Holstein-Primakoff transformation: 



a k — v k b k — fi k b*_ k 
a -k = Vkb*_ k - \x h b k , 



(23) 



where 



v k = cosh 



Xk 
2 ' 



/i fc = e i ' 3fe sinh^, 
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C k + ReCfc 
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1/2' 



Substituting Eq. (23) into Eq. (20), we obtain 
d t b k - \{u k + irj k )b k = 0, 



(24) 



(25) 



where 



u 2 k = (C k + ReC k ) 2 - (1 + A 2 ) \D k e 2i ^ +C k 
Vk = A(C fe + Re(7 fe ) 



(26) 

(27) 



Equations (26) and (27) express the dispersion relation and a damping rate, 
respectively. 



3 The instability threshold 



Now we consider the case where the microwave field h cos ut is applied. The 
equation of motion of a k corresponding to Eq. (20) is 

d t a k + i(l — \\)A k a k + i(l — i\)B k a*_ k — i(l — iA)^ cos cu p t ■ a k = 0. (28) 



After the the Holstein-Primakoff transformation (23), we substitute the fol- 
lowing equations into Eq. (28), 



bh(t) = b° h (t)exp[i(w p /2)t-7} h t], 
b*- k (t) = b°_* k (t) exp[-i(wp/2)t - Vk t], 



(29) 



since the resonance occurs at uo k = uj p /2. Then the slowly-varying variable b k 
satisfies 



d\b 



2uo 
k 



UJ, 



u k 



\pk\ 



K = o, 



(30) 



where 



\p h \ = uj h Vl + A 2 



D k e 2i ^+C k 

2uu 



(31) 



Therefore, the exponentially increasing solution for b k is 

b k oc exp [flpfcl -rj k )t + i(u> p /2)t] , (32) 

where 

IPfel > Vk- (33) 

The instability threshold a;" 1 * is now given as 

<" = -A, mini - J* -, I . (31) 



VI + A 2 * D k e 2i ^+C k 



The threshold, Eq. (34), is obtained by using only the LL equation, but cannot 
explain the experimentally-observed instability threshold for parallel pumping. 
In fact, Eq. (34) proves to give an instability curve totally different from the 
real butterfly curve. Since the relaxation phenomenon is essentially nonlinear 
one, the linearizion analysis is not sufficient to discuss the instability threshold. 
The problem is beyond the purpose of this paper, and it will be discussed in a 
forthcoming paper [9]. This difficulty, however, can be overcome by replacing 
the damping rate r\ k by a suitable spin-wave line width. Using the spin-wave 
line width AH k , we rewrite Eq. (34) as 



o;? rit = u n min 



AH k 



p 



k 



D k e^"C k 



(35) 



Here we adopt a simple trial AH k function [3], 

AH k = A + A 1 sm 2 (26 k ) + A 2 k, (36) 

where Aq, A\ and A 2 are adjustable parameters. 
4 Butterfly curves 



Typical butterfly curves of the threshold o;™ 1 * have a cusp at a certain static 
field: as the static field increasing, uj^ lt decreases below the cusp point and 
increases above that. For static fields below the cusp point, the minimum 
threshold modes corresponds to a spin wave propagating with 6 k = ir/2. As 
the static field increases, the wave vector k of the threshold modes decreases, 
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Fig. 2. Theoretical butterfly curves for d 
angle a. 
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and k ~ at the cusp. For static fields above the cusp point, the wave number 
remains at k ~ and 6^ decreases from tt/2 to 0. 

Figure 2 shows theoretical butterfly curves calculated by using Eqs. (35) and 
(36). The material parameters used in the calculation are typical values for 
yttrium iron garnet (YIG) materials: J7| = 1.77 x 10 7 rad/(s-Oe), 47rMo = 
1.75 x 10 3 Oe, D = 5.4 x 10 _9 Oe-cm 2 /rad 2 . The other parameters are as follows: 
uj/(2tv) = 9.5 x 10 9 Hz, d = 15^m, A = 3.0 x 10- 2 /(4irM ), AJAq = 0.4, 
A2M0 = 3.0 x 10~ 4 /d We calculate the static field value for the cusp point 
with k = and 9j, in a bulk approximation: we assume fj~ = 1 in Eq. (22). 
Setting Uk = ou p /2 in Eq. (26), we obtain the static field for the cusp point 



to 



.crit 
H 



cos 2 a 



1 1 

2 + 2 



u* + l. 



(37) 



While the cusp point shifts to lower static fields as the oblique angle a in- 
creases, the typical shape of curves is found to show little noteworthy change. 

In some cases, however, curious butterfly curves can emerge. Figure 3 is a 
theoretical butterfly curve for the case of d = 5/im and a = n/2. The other 
parameters used in the calculation are the same as used in Fig. 2. There appear 
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Fig. 3. Theoretical butterfly curve for d = 5/um. The external field is applied parallel 
to the film plane. 

multiple cusps in the low field region of the butterfly curve. These cusps are due 
to different standing spin-wave modes across the film thickness. Let us recall 
that the z component of wave vectors are quantized. Effects of quantization 
are essential when the boundary conditions become important. Comparing the 
results in Figs. 2 and 3, one may say multi-cusp feature of a butterfly curve can 
be seen when the thickness of a film is very small. This finding is supported 
by some experimental studies [4,5]. In Ref. [5], a standard butterfly curve of 
the instability was shown for 15.9-//m-thick yttrium iron garnet (YIG) film 
under in-plane external field. However, in Ref. [4], multiple butterfly curves 
appeared for a YIG film under in-plane external field. The difference between 
the two is just the thickness. The thickness of the film in Ref. [4] was 0.5/zm. 



5 Discussion 



We have revealed how butterfly curves depend on the oblique angle between 
the external field and the film plane. We have calculated theoretical butterfly 
curves by using the LL equation together with the spin- wave line width AH},. 
The calculation was performed under the assumption that parameters Aq, 
A\ and A 2 for AHk do not depend on the oblique angle. These parameters 
were originally introduced to fit a theoretical butterfly curve to experimental 
data [3], and might depend on the oblique angle. 

We have also shown qualitative features of the novel aspect of butterfly curves 
with multiple cusps for parallel pumping. Those multi-cusp curves come from 
quantized standing-wave modes across the film. However, it is not clear how 
such cusps appear, since there remains an ambiguity about how to evaluate 
A Q , Ai and A 2 . Figure 3 is one of possible theoretical curves. To proceed to 
quantitative evaluations, the experimental test is highly desirable to confirm 
the features predicted here. 
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A Demagnetization field 



First, we should use the (x, y, z) coordinate system rather than (x, y, z) to use 
Eqs. (15) and (16) in the Poisson equation (4), and employ the transformation: 
V 2 — > V 2 , d x — ► cosadx — sinac^. We consider only the case of — \ < z < \. 
For even modes, k z = 2mix (m: integer), the Poisson equation is rewritten as 



V 2 $ = ^e 

k 



i(k x x+k y y) 



[ik- cos k z z + k z sin a sin k z z)a] c 



+ (ik+ cos k z z + k z sin a sin k z z)a*_ 



(A.1) 



where k + = k x cosa + \k y and A;_ = ^ cos a — \k y . For odd modes, k z 
(2m + 1)tt (m: integer), 



V 2 $ = ]Te 
k 



i(k x x+k y y) 



(ik- sin k z z — k z sin a cos k z z)ak 
iA;_i_ sin k 7 z — k 7 sin a cos k?z)a* u 



(A.2) 



Here we note a well known fact: the Poisson equation, 



V 2 (/>(r) = -Anp(r) 
has the solution as 



(A.3) 



p(r') dV 



\r — r 



(A.4) 



Applying Eq. (A.4), we obtain the solution of Eqs. (A.l) and (A.2): for even 
modes, 
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+ k z sina ■ (oi, + a 



dV 



fe l -fc^ / if_f| 



3 i(fes'+M') sin j^ 1 ; ( A5 ) 



for odd modes, 






dV' 



+ "-*•/ / \~t 



Ak x x>+k y §') gin r ~, 



dV' 



- k z sin a • (Ofe - a* r ) / —, e^*'*^ cos k z ~z' ) (A.6) 



°fc w -fc^ / |f.'_ r \ 



After integration, Eqs. (A. 5) and (A.6) become 



g(k x x+k y y) 

Hr) = ~ E p {i(fc-a* + MI*) 



cos Ls - (-l) m e~^ /2 cosh fc. 



+ A; 2 sin ct • (o^, + a%' 



sin ~k z ~z - (-l) m e"^ /2 b- sinh fc ± 5 

fci 



(A.7) 



and 



e i(k x x+k y y) _ _ . _ 

$(f) = - E ^ {i^-Ofc - MI*) [sin k z z - (-l) m e"^ /2 sinh k ± z 



k z sin a ■ (a^ — a*? ) 



cos M + (-l) m e-^/ 2 ^ cosh A; ± 5 



,(A.S 



respectively. 



B Equation of motion for a^ 



First, let us calculate the derivative of $ in Eq. (19). For even modes, k z = 2mn 
(to: integer), 



(cos adx + \dy — sin ad z )<& = E 



j(k x x+k y y) 
P 



X 



{-(fc+fe-Ofc + fe+o^) cosk z z- (-l) m e fex/2 cosh /c_lz 

+ i^ + ^ sin a ■ (o^ + a*_^) sink z z — (—l) m e~ k± ^ 2 (k z /k±)smh.k_i_z \ 
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+ sin a | i{k_a k + k + a*_ 



+ u ~y 



-k sin k~z - (-l) m e~' k ^ 2 k , sinh k 



+ k 2 z sina- (a k + a*_ k ) cos k z z - (-l) m e fc±/2 cosh k±z } 



(b.i; 



where A; + 

(2m + l)7r (m: integer 



k x cos a + ifcy and k. 



k x cosa — iky. For odd modes, k z 



(cos ad x + idy — sin ad z )<& = ^ 






X 



{-(hk-a~ k - ~k 2 + a*_~ k ) 



ik+k z sin a ■ (a k + a*_ 
at — k+a* 
+ k z sin a • (a k — a*_ k ) 



sin k z z - (-l) m e k±/2 sinh A;_l,2 

tos k z z - (-l) m e-~ k±/2 (k z /k±) cosh k ± z 



+ sin a | i(^-Wfe — *-+«_£, 



k z cosk z z - (-l) m e k±/2 k± cosh k±z 
sin A; 2 5 - (-l) m e~ fe±/2 sinh A;_l5 } 



(B.2) 



Let us project Eqs. (B.I) and (B.2) onto cos k z z and sin k z z, respectively. For 
—1/2 < z < 1/2, the projection F(z) of a function f(z) onto cosA^i: is given 
by, 



1/2 1/2 

F(z) — / dzf(z) cos k z z / / d5cos 2 A; z 5. 

-1/2 ' -1/2 



(B.3) 



After the projection, Eqs. (B.I) and (B.2) are reduced to 



(cos ad x + idy — sin ad z )<& = y^ 

k 



j(k x x+k y y) 



1 - (2 - 6, 



mOj 



l-e 



where &,■ is the Kronecker delta, and 



A; 2 1 \ ^ A; 2 

x •{ {k + k- + A; 2 sin 2 a)^ + (&;+ + A; 2 sin 2 a)a*_ k \ cos k z z, (B.4) 



e i(fe I s+fc !/ j/) r 

(cosa<9 £ + i(9y — sinadz)$ = ^ n 1 1 — ^ (l 

k { 



k 2 



x 



I (k + k- + k 2 sin 2 a)a k — (&;+ + A; 2 sin 2 a)a*_ k \ sin A; z z, (B.5) 



respectively. Let us impose a restriction on the kinds of paramete rs about 
wave vectors to describe the equation of motion: we only use k, k± = Jk 2 + A; 2 

k z = nn {n: integer), 6 k and ip k - Here, sin# fc = A;_|_/A;; cosy?*, = k x /k±; sinip k = 
k y /k±. The restriction leads to the following: 
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k + k_ + k z sin 2 a = k\ + 2(k± cos (p k + fc 2 sin a) k z sin a; 
£;+ + k 2 z sin 2 a = A;^ e 2i ^ +2(A; ± e 1Lp k + k z sin a)A; z sin a; 



k ± = ^Jk 2 -k 2 z . (B.6) 

From the above two equations and Eq. (19), we obtain Eq. (20): 

d t a k + i(l - iA)A fc a fc + i(l - i\)B k a*_ k = 0. (B.7) 

For even modes, 

A k = C k + ReC fe , B k = (D k e 2i ^ +C k ) , (B.8) 



and for odd modes, 

A k = C k + ReC k , B k = - (D k e 2i ^ +C k ) . (B.9) 

Here we note that C k , C k and D k are described in Eq. (22). The difference 
between Eqs. (B.8) and (B.9) is just the sign of B k . Therefore, for both even 
and odd modes, we may write 

B k = (-ir(D k e 2 ^+C k ). (B.10) 

References 

[1] E. Schlomann, J. J. Green and U. Milano, J. Appl. Phys. 31, 386S (1960). 

[2] H. Suhl, J. Phys. Chem. Solids 1, 209 (1957). 

[3] See, for example, M. Chen and C. E. Patton, in Nonlinear Phenomena 
and Chaos in Magnetic Materials, edited by P. E. Wigen (World Scientific, 
Singapore, 1994), pp. 33-82. 

[4] B. A. Kalinikos, N. G. Kovshikov and N. V. Kozhus, Sov. Phys. Solid State 27, 
1681 (1986). [Fiz. Tverd. Tela (Leningrad) 27, 2794 (1985).] 

[5] G. Wiese, L. Buxmzn, P. Kabos and C. E. Patton, J.Appl. Phys. 75, 1041 
(1994). 

[6] P. Kabos, M Mendik, G. Wiese and C. E. Patton, Phys. Rev. B 55, 11457 
(1997). 

[7] M. Lakshmanan and K. Nakamura, Phys. Rev. Lett. 53, 2497 (1984). 

[8] F. J. Elmer, Phys. Rev. B 53, 14323 (1996). 

[9] K. Kudo and K. Nakamura, in preparation. 



14 



